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ABSTRACT 

The  consistency  and  asymptotic  normality  of  a  linear  least  squares 
estimate  of  the  form  (X*X)~X*Y  when  the  mean  is  not  X£  is  investigated  in 
this  paper.  The  least  squares  estimate  is  a  consistent  estimate  of  the 
best  linear  approximation  of  the  true  mean  function  for  the  design  chosen. 
The  asymptotic  normality  of  the  least  squares  estimate  depends  on  the  design 
and  the  asymptotic  mean  nay  not  be  the  best  linear  approximation  of  the  true 
mean  function. 

V 

Key  words:  Consistency,  asymototic  normality,  best  linear  approximation, 
model  robustness. 
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Linear  Least  Squares  Estimates 
and  Nonlinear  Means 

1.  Introduction 

In  the  standard  linear  models  theory,  the  mean  of  an  observation  vector 
Y  is  given  Xg  where  X  is  a  known  matrix  of  constants  and  g  is  an  unknown 
vector  of  parameters.  The  estimate  8  =•  (X'X)  X'Y,  where  (X'X)  is  a  gen¬ 
eralized  inverse  of  X"X,  has  many  well-known  and  desirable  properties.  Con¬ 
sidering  that  both  regression  and  analysis  of  variance  problems  are  part 
of  the  linear  models  theory,  estimates  like  8  are  among  the  most  widely  used 
of  all  parameter  estimates. 

A 

If  the  nean  of  Y  is  not  X8,  the  properties  of  8  are  not  described  by 
the  standard  linear  models  theory.  Indeed,  it  may  not  be  clear  what  (if 
anything) is  estimating.  The  behavior  of  8  for  large  samples  when  the  mean 
is  not  XB  is  investigated  in  this  paper.  This  behavior  is  found  to  depend 

A 

on  the  design  used  in  the  experiment.  But,  for  a  given  design,  8  is  shown 
to  estimate  the  best  linear  approximation  of  the  true  nean  function  in  a  sense 
to  be  defined. 

In  Section  2,  the  model  and  notation  are  defined.  In  Section  3,8  is 
shown  to  be  a  consistent  estimate  of  the  best  linear  approximation.  The 
asymptotic  normality  of  8  is  investigated  in  Section  A.  Section  5  contains 
some  examples. 

Other  authors  have  considered  questions  related  to  those  addressed  herein. 
The  question  of  model  robustness  what  8  is  estimating,  if  the  mean  is  not 
XB,  has  been  investigated  by  authors  such  as  Box  and  Draper  (1950),  Atwood 
(1971),  Stigler  (1°71)  and  Huber  (l',75).  These  authors  considered  specific 
alternative  mean  functions  to  X8.  Royall  and  Herson  (1973)  have  described 
the  mean  of  8  for  an  arbitrary  mean  function.  Surprisingly  recently  other 
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authors  such  as  Drygas  (1976),  Lai,  Robbins  and  Wei  (1978),  Anderson  and 
Taylor  (1979)  and  Wu  (1989)  have  considered  the  consistency  of  8  under 
various  conditions,  but  always  under  the  assumption  that  the  mean  is  X(3. 


2.  ftodel  and  Notation 

Let  X  denote  a  subset  of  an  r-dinensional  Euclidean  space  (r  a  l). 

X  is  the  set  of  possible  values  of  the  independent  variable  x.  A  design  for 

a  sample  of  size  n  is  a  specification  of  n  points,  x,  ,  from  X 

l  ,n  ‘  n,n 

where  the  n  points  specify  the  values  of  the  independent  variable  at  which 

observations  are  to  be  taken.  The  x.  need  not  be  all  distinct.  If  a 

i,n 

point  x  is  repeated  p  times  then  p  observations  are  taken  at  x.  A  design 

can  be  completely  described  by  a  discrete  probability  measure  on  X 

where  the  probability  assigns  to  a  point  x  is  the  proportion  of  the  n 

observations  to  be  taken  at  x.  5  will  also  be  called  the  design. 

n 

The  observation  vector  is  Y  *  (Y,  ,...,  Y  )'.'.  It  is  assumed  that 

n  l,n  n,n 

for  each  n 


i.n 


=  m(x.  _  )  +  e. 


i.n 


i,n 


i  =  1, . . . ,n. 


where  n,  the  unknown  mean  function ,  is  a  real-valued  bounded  and  continuous 

function  on  X  and  e.  ,  n  *  1,2,...;  i  =  l,...,n,  are  identically  distributed 
l  ,n 

2 

random  variables  with  mean  zero  and  variance  o  .  It  is  also  assumed  that 


e,  ,,..,e _  are  independent  for  each  n  »  1,2,..,  .  Thus  Y,  ,...,Y 

l,n  n,n  i,n  n,n 

2 

are  indenendent  random  variables  with  13  Y.  =  m(x.  )  and  variance  o  . 

i,n  i,n 

let  f(x)  «  (fj(x),...,  f  (x))'  denote  a  p  *  1  vector  of  real-valued, 

bounded  and  continuous  functions  on  X.  For  a  given  design,  x.  ,...,x  , 

i ,n  n>n 

let  be  the  n  *  p  matrix  with  (i,  j)th  element  f^(x^  n).  The  asymptotic 

behavior  of  the  estimate  *  (X'X  )~X'Y  will  be  investigated  in  this  »aper. 

n  n  n  n  n 


8n  will  be  called  a  linear  least  squares  estimate  since,  for  a  given 

observation  y  *  (y,  ,...,y  )',  8  is  the  vector  8  which  minimizes 

n  l,n  n,n  n 

E?  .  (y.  -  B'f(x.  ))2.  Conditions  under  which  X'X  is  non-singular  and 

1=1  i,n  i,n  n  n 

is  uniquely  defined  are  given  in  Lemmas  3.1  and  3.4. 

Let  C  be  a  probability  measure  on  X.  A  p  *  1  vector  B(n,£)  will  be  called 
a  best  linear  approximation  of  mfx)  if 

/ (m(x)  -  B'(m,£)f(x))2d5(x)  =  inf  /( n(x)  -  8'f(x))2d£ (x)  .  (2.1) 

8 

Since  8(m,S)  depends  on  the  unknown  mean  function  m,  B(n,£)  is  a  parameter. 

Let  M(C)  denote  the  p  x  p  matrix  with  (i,j)th  element  Jf^(x)  f  ^  (x)  d£(x)  and 
let  c(£)  denote  the  p  x  1  vector  with  (i)th  coordinate  /f^(x)m(x)  d£(x).  In 
optimal  design  theory  (see,  e.g.,  Kiefer  (1962))  a  multiple  of  M(£)  is  called 
the  information  matrix  of  £.  If  H(Cn)  is  non-singular,  a" n  M  (£)  is  the 

A 

covariance  matrix  of  By  equating  the  partial  derivatives 

a/(m(x)  -  B'f(x))“d£(x)/36j ,  j  =  l,...,p,  to  zero  it  is  easily  verified  that 

if  M(g)  is  non-singular  then  8(m,£)  is  unique  and  equals  M  *(£)c(£).  Under 

a 

conditions  relating  a  sequence  of  designs  to  £,  it  will  be  shown  that 
is  a  consistent  estimate  of  B(m,£)  and  6n  is  asymptotically  normal  with  mean 
B(m,£).  In  this  sense,  the  best  linear  approximation  8(m,£)  is  the  parameter 
being  estimated  by  the  linear  least  squares  estimate  Bn  when  the  mean  m(x)  is 
not  necessarily  linear. 

3.  Consistency 

Let  £n  be  a  fixed  sequence  of  designs  for  sample  sizes  n  =  1,2 . 

Throughout  it  will  be  assumed  that  converges  weakly  to  a  probability  measure 


C  where  weak  convergence  is  defined  in  Billingsley  (1968).  The  main  result 
in  this  section,  the  proof  of  which  is  deferred  until  the  end  of  the  section, 
is  the  following  theorem. 

Theorem  3.1:  If  ?i  (£)  is  non-singular,  then  is  a  consistent  estimate 
of  8(m,£)  in  that  the  random  vectors  sn  converge  in  probability  to  the  real¬ 
valued  vector  8(m,£)  =  fi”*(£)  c(£). 

It  is  important  to  note  that  the  parameter  £S(m,£)  which  Bn  estimates 
depends  on  £.  The  experimenter  chooses  £  when  the  sequence  of  designs 
£j,  *s  ch°sen*  5  should  be  chosen  so  that  the  definition  of  best 

linear  approximation  in  (2.1)  accurately  reflects  how  the  experimenter  wishes 
to  measure  the  closeness  of  g'f(x)  to  m(x).  If  the  support  of  £  consists 
only  of  a  finite  number  of  points,  g"f(x)  will  be  compared  to  m(x)  only  at 
these  points  in  determining  g(n,£).  If  the  aim  is  to  estimate  a  B  so  that 
B'f(x)  is  close  to  m(x)  for  all  x  in  X,  a  choice  of  a  £  whose  support  is  all 
of  X,  e.g.,  uniform  on  X,  seems  more  appropriate.  In  optimal  design  theory 
(see,  e.g.,  Kiefer  ( 1C 50)  or  Karlin  and  Studden  (1966))  the  optimal  design 
often  has  a  support  with  a  finite  number  of  points.  These  designs  may  not 
be  very  appropriate  if  the  mean  m(x)  is  not  of  the  form  B'f(x).  The  fact 

A 

that  what  B  is  estimating  depends  on  the  design  if  m(x)  i  B'f(x)  has 
n 

been  recognized  previously.  See,  for  example.  Draper  and  Smith  (1966,  Chapter 
2,  Section  12)  or  Royall  and  Herson  (1973). 

The  following  lemma  gives  a  condition  under  which  M(£)  will  be  non¬ 
singular.  The  functions  fj (••),...  ,f  (x)  are  called  linearly  independent  -  £ 
if,  for  a  p  x  1  real-valued  vector  a,a'f(x)  *  0  a.s.  £  imnlies  a  *  9. 

Lemma  3. 1 :  If  f^x) ,. . .  ,f^(x)  are  linearly  independent  -£  then  M(£) 


is  non-singular. 
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Proof:  Let  a  be  a  p  *  1  real-valued  vector.  Assume  M(^)a  =  0*  To 
show  11(C)  is  non-singular  it  suffices  to  show  a  =  0.  M(£)a  *  0  implies 

0  *  (£)<*  *  /(a'f(x))2dS(x) .  This  implies  a'f(x)  =  3  a.s.  £ .  Since 

fjCx) ,  f  (x)  are  linearly  independent  -£,  a  =  0.  j| 

To  study  the  asymptotic  behavior  of  it  is  useful  to  note  that 


6  *  (X'X  )  X'Y 

n  n  n  n  n 


(3.1) 


where  e  *  (z  e  „)'  and  n  *  (n(x.  „),...,  m(x  ))'  .  The  asymptotic 

n  1  ,n  n,n  n  l,n  n,n 

behavior  of  n  *X'e  .  n  *X;*m  ar.d  n  ^X'X  is  described  in  the  following  three 
n  n  n  n  n  n 

lemmas . 

Lemma  3.2:  Under  the  model  the  random  vectors  n-1X'en  converge  in 
probability  to  0. 


Proof:  It  suffices  to  show  that  for  any  6  >  3, 


lim  P(|z"b1  fr(xi  n)Ci  |  >  n  j)  =  0  for  r  =  l,...,p. 


(3.2) 


2 

Let  a  =  max  sup  f  (x) .  Since  the  f  are  bounded,  a  <  <».  Fix 
isr<p  y  T  r 

5  >  0.  By  Chebyshev's  Inequality  and  the  fact  that  e,  , . . .  ,e  are  uncorrelated 

1  jii  n,n 


with  mean  zero. 


P(h".,  fr(»lin)«ifn)2/n2j2 

°  ^*1  ^r^i.n^0"5  *  <?2a/n62.  Thus  Equation  (3.2)  is  true.  1 1 

For  the  sake  of  completeness  it  should  be  noted  that  the  result  of  Lemma 

3.2  and  hence  the  consistency  result  of  Theorem  3.1  holds  under  the  condition 

that  for  each  n,  e1>n> • • • are  uncorrelated  with  means  all  zero  and  variances 

all  bounded  above  by  o2.  The  stronger  condition  that  e.  ,...,e  are  i  i  d 

i »n  n,n  ’  ’  * 

is  not  necessary. 


Lemma  3.3:  Under  the  assumptions  of  the  model,  lim  n  *X'm^  =  c(£). 
Proof:  Coordinatewise,  this  result  is 

lim  n'1  I1}  f  (x  )m(x.  )  =  /  f  (x)ra(x)d£ (x) ,  r  =  l,...,p. 

n-*»  i=l  r  i,n  i,n  r 

Note  that  n"*  £?  ,  f  (x.  )m(x.  )  =  /f  (x)m(x)d£  (x) . 

1=1  rv  i,n  i,n  1  r  n 

Since  f  and  m  are  bounded  and  continuous  and  £  converges  weakly  to  £ 
r  n 

(3.3)  is  true. | | 


(3.3) 


Equation 


-1. 


■'*(£)  * 


Lemma  3.4:  If  f?(£)  is  non-singular  then  (i)  lim  n  X'X 

_  _ _  n  n 

n-**> 

(ii)  n_1X"X  is  non-sinpular  for  all  sufficiently  large  n  and 

n  n 

(iii)  lim  (n'Vx  )“  =  ifJ(£).' 

11 

Proof :  Elementwise,  statement  (i)  is 

lim  n'1  Zn  .  f  (x.  )  f  (x.  )  =  / f  (x)f  (x)dC(x),  r,  s  =  l,...,p.  (3.4) 

1=1  r  1>n  s  1."  r  s 

Note  that  n-1  E?  f  (x.  )f  (x.  )  =  / f  (x)f  (x)d£  (x) . 

i=l  r  i,n  s  i,n  r  s  n 

Since  f  and  f  are  bounded  and  continuous  and  £  converges  weakly  to  €, 
r  s  n 

Equation  (3.4)  is  true. 

Since  the  determinant  is  a  continuous  function,  by  (i), 

lim  |n_1X'X  |  =|  M(£)|  ^  0.  Thus  (ii)  is  true. 
n  n  1 

n>  _1  _1 
If  A  and  A  are  non-singular  matrices  and  lim  A  =  A,  then  lim  A  =  A  . 
n  6  n  n 

j  r.-*»  n-*» 

By  (i)  and  (ii) ,  n  X'X  is  non-singular  for  all  large  n  and  (iii)  is  true.|| 

n  n 

Proof  of  Theorem  3.1:  This  theorem  follows  from  Equation  (3.1)  and  Lemmas 
3.2,  3.3  and  3.4. | j 
4.  Asymptotic  Normality 

1/2  ' 

In  this  section,  the  asymptotic  normality  of  n  (g^  -  0(m,£))  is  investi¬ 
gated.  As  in  Section  3,  it  is  assumed  that  £n  is  a  fixed  sequence  of  designs 
for  sample  sizes  n  =  1,  2,...  whi«h  converges  weakly  to  a  probability  measure 


£  .  The  asymptotic  normality  result  is  given  in  Theorem  4.1. 


s 

Theorem  4,1:  Assume  M(£)  is  non-singular.  Let  m  =  (m(x,  m(x  ))r 

~ ~ n  l,n  n,n 

1/2  — 

Assume  that  lira  n  ((X'X  )  X"n  -  B(m.,f;))  exists  and  equals  the  p  x  1  vector 
n-*»  nr*  n  n 

1/2  " 

b.  Then  n  (f$  -e(m,^))  converges  weakly  to  a  multivariate  normal  random  vector 

with  mean  b  and  covariance  matrix  o  K  (C). 

Proof:  Let  T  =  n  ^2X'e  .  Note  that 
-  n  n  "> 

1/2  *  -1  —  1/2  — 
n  (*n  -  B(».C))  -  (n  AX^Xn)  T^  „  '"((X^)  X^  -  *(nfO).  (4.1) 

By  assumption,  the  last  term  converges  to  b.  By  Lemma  3.4,  lim  (n  *X'X  )  =  I!  *(£). 

n-*°°  n  n 

So  to  prove  the  desired  result  it  suffices,  by  an  extension  of  Slutzky’s  Theorem 

(Billingsley  (1968),  Theorem  5.5),  to  show  that  the  random  vectors  Tn  converge 

weakly  to  U  where  U  is  a  p-dimensional  multivariate  normal  random  vector  with 

2 

mean  0  and  covariance  matrix  o  f'(fj. 


To  prove  the  convergence  of  T  to  !j,  it  suffices,  by  the  Cramer-Wold  device 

n 

(Billingsley  (1962),  p.  4G) ,  to  show  that  ot'T^  converges  weakly  to  a'U  for  all 

- 1/2 

p-dimensional  vectors  a.  Fix  a  /  0.  Let  =  n  (a'f(x.  ))e.  •  Then 

3  ,n  3  >n  3  »n 

a'T  =  rn  ,  so  it  suffices  tc  shew  that  I1?  ,  W.  converges  weakly  to  a'U. 

n  j-1  ]  ,r.  3  =  1  1  ,n 

are  independent  random  variables  with 


n,n 


By  the  definition  of  t.  ,  h 
'  J,»*  l,n 

.  i  ?  2 

zero  means  and  variances  respectively  equal  to  r.  *(ct'f(x..  ^))-a“.  Thus  it 


suffices,  by  the  Normal  Central  Limit  Theorem  (Lo6ve  (1963),  p.  288),  to  show  that 

lim  In  ,  LV'2  =  aV':(C)a  (4.2) 

n->«  j  =  l  j  ,n 

and  that 

lim  E1?  ,  E ’ .  I(| W.  |>  5))2  =  0  for  all  5  >  9  (4.3) 

3  =  1  j  ,n  1  3 ,n'  ' 

where  I  is  the  indicator  function. 

To  verify  Equation  (4.2)  note  that  E°  ,  EW?  =  a “a' (n  *X'X  )  a .  By  Lemma 
1  1  j  =  l  3  ,n  n  n 

3.4,  lim  (n  *X'X  )  =  >  1(0  so  Equation  4.2  is  true. 
ii'Too  n  n 

Finally  to  verify  Equation  4.3,  let  E  =  synja'f(x)|.  Since  the  f^  are 

bounded,  B  <  » .  Fix6  >  0.  Since e  ,  , . . .  ,e  are  i.i.d.  random  variables , 

l,n  n,n 


0 


ew2  i  ( 

IW.  1  >  6 

) 

1  .n 

1  3  ,n 

lim 

-1  _n 
n  l . 

,Ee?  I( 

1  e  • 

3s 

'1  3  >n 

1  3 

lim 

E*?.i 

'‘I'i.iI 

> 

55  °*ll 


In  theorem  4.1,  the  asymptotic  mean  b  depends  on  the  unknown  function 

m(x).  In  order  for  Theoren  1.1  to  be  useful  for  making  large  sample  inferences 

about  say  to  construct  confidence  sets  for  e(m,f;),  b  must  be  zero  for 

all  bounded,  continuous  functions  ir.(x).  Lemmas  3.3  and  3.4  show  that  (X'X  )  x'm 

n  n  n  n 

converges  to  f$(m,£)  for  all  bounded,  continuous  n(x).  It  is  reasonable  to  expect 

that  b  =  0  if  £  converges  to  S  fast  enough  in  some  sense.  Conditions  under 

which  this  is  true  are  examined  in  the  remainder  of  this  section.  A  useful 

result  in  this  regard  is  given  in  Corrollary  4.1.  It  should  be  noted  that  one 

important  situation  in  which  b  =  0  is  if  m(x)  =&'f(x)  for  some  6  .  In  this 

case  (X'X  )~ X'm  =  6  =  6(n,Q  for  all  large  n. 

n  n  n  n 

To  simplify  notation,  for  the  remainder  of  this  section  assume  that  X  =  (u,v) 
an  interval  on  the  real  line.  Similar  results  hold  if  x  is  a  r-dimensional 
vector.  Let  £^(x)  and  £(x)  also  denote  the  distribution  functions  of  the  pro¬ 
bability  measures  ^  and  C.  Assume  there  is  a  function  h(x)  such  that 
1/2 

lim  n  (£^(x)  “  £lx))  =  almost  everywhere  with  respect  to  Lebesque 
measure  on  (u,v).  Assume  that  sun  r.1  "|Cn(x)  -  £(x)jdx  <  00 .  Assume  the 
functions  m,  f  , . . .  ,  f  are  differentiable  on  (u,v).  Let  g'  denote  the  deriva¬ 
tive  of  a  function  g.  The  mean  vector  b  in  Theoren  3.1  can  be  written  in  terms 
of  the  function  h.  To  obtain  this  result,  the  following  lemma  will  be  used. 

Lemma  4.1:  Let  <}>(x)  be  a  real  valued  function  defined  on  (u,  v)  with  a 
derivative  ^'and  let  f.(x)  be  a  distribution  function  with  support  contained  in 


(4.4) 


$(x)dG(x)  =  <j>'(x)  (1  -  G(x))dx  +  $(u). 

Proof:  (j>(x)dG(x)  =  f ^  /*<♦)- (u)du  dG(x)  +  <*>  (u) . 

Equation  4.4  follows  from  Futini  s  Theorem  (Lo6ve  (1963),  p.  135)  by  inter¬ 
changing  the  integration  order. j j 

The  next  theorem  gives  an  expression  for  the  near,  b  of  Theorem  4.1  in  terms 

of  matrices  which  depend  on  the  function  h.  Let  d  denote  the  p  *  1  vector  with 
fV 

(i)th  coordinate  -J  u  (f^  (x)rn(x))  (x)dx.  Let  P  denote  the  p  *  p  matrix  with 

(i,  j)th  olorent  /^.  (x)f.(x))'i(x)dx. 

Theorem  4.2:  Assume  ."(£}  is  non-singular .  Then  b  =  fi  *(C)(d  +  DM  *(£)c(C)). 

Proof:  Let  A  =  n  *X'X  .  "y  Lemma  3.4,  A  is  non-singular  for  all  large 
n  r*  n  n 

n  so  we  shall  write  A  1  for  A- 

n  n 

\  fj  _  1  _  1  _  1 

b  =  lim  n  (Anx(n  -  M  ‘(Oc(«) 

=  Ai22  nl/2(An1(n'lxA  •  c(5))  +  {\!  -  v_1(5))c(C)). 

To  prove  the  desired  result  it  suffices  to  prove  that 

lim  n1/2A~X (n-1X'm^  -  c(C))  =  M_1(C)d.  (4.5) 

and 

lim  °1/2  (A’1-  f!_1(5))c(C)  =  M'1(C)Df.;'1(C)c(C).  (4.6) 

By  Lemma  3.4,  lim  A  *  *  M  *(£).  Thus  to  prove  Equation  4.5  it  suffices  to 
n-H»  n 

show  that 

lim  nl/2(/'^  fr(x)m(x)dCn(x)  -  J*  fr(x)m(x)dC(x)) 

■  -lim  fVu  (fr(x)m(x))'n1/2(Cn(x)  -  C(x))dx  (4.7) 

~  -  I u  (fr(x)m(x)) 'h(x)dx,  r  =  l,...,p. 

The  first  equality  is  true  by  Lemma  4.1  and  the  second  equality  is  true  by  the 

1/? 

assumptions  made  about  n  "(Cn(x)  -  c(x))  an<l  the  Dominated  Convergence  Theorem 


(Loeve  (1963),  p.  125). 

Note  that  A'1-  h'1^)  =  A*1  (  (C)  -  A^f-f1^).  3y  Lemma  3.4, 
lir  A^  =  *(£)•  Thus  to  prove  Equation  4.6  it  suffices  to  show  that 

1,1/2  (/u  fr(x^fs(x)  d5(x)  '  /u  fr(x)fs(x)  d?n(x)) 

=  lyn  /JJ  (fr(x)fs(x))"nx/-  (Cn(x)  -C(x))dx  (4.8) 

*  I*  (f  (x)f.(x))h(x)dx,  r  =  1, . . . ,p;  s  =  1, .  . .  ,  o. 

u  r  o 

The  first  equality  is  true  by  Lemma  r.i  and  the  second  equality  is  true  by  the 

1/2 

assumntions  made  about  n  (C  (x)  -  C(x))  and  the  Dominated  Convergence  Theorem.)  | 

1/2 

Corrollary  4.1:  Assume  life)  is  non-singular.  If  lim  n  (C  (x)  -  E(x))  =  0 
— -  n-Ko  n 

1/2  * 

almost  everywhere  with  resuect  to  Lebesque  measure  on  (u,  v)  then  n  ' (p  -  g(m,c)) 

converges  weakly  to  a  multivariate  normal  random  vector  with  mean  0  and  covariance 
2,  -1 , 

matrix  a  , .  (C) . 

i  /a 

Proof :  h(x)  =  lir.  n  7~(Cn(x)  -  C(x))  =  0  a.e.  implies  d  =  0  and  D  =  0. 

Thus  by  Theorem  4.2,  the  mean  b  in  Theorem  4.1  is  O.j) 

5.  Examples 

In  this  section  two  examples  are  considered.  In  the  first  example  a  fairly 
general  method  of  constructing  designs  which  satisfy  the  conditions  of  Corollary 
4.1  is  given.  The  second  example  gives  a  sequence  of  designs  for  which  the 
function  h(x)  is  non-zero. 


Example  1 :  Let  c(x)  he  a  fixed  continuous  and  strictly  increasing  distri¬ 
bution  function  on  the  real  line.  Let  C  * (x)  denote  the  inverse  of  C(x) .  Let 


n^,  n.,,...  and  n^t  r.., , . . .  be  sequences  of  positive  integers  such  that  n^/nn 

l/7 

is  an  integer  for  each  i  and  lim  r../n.  ~  =  0.  For  k  =  1,...,  n./m.  ,  let 

i-*«=  i  i  ii 

I,  ^  denote  the  interval  (E  *((k  -  l)rn/m),  E  *  (km Am  ) ) .  Let  Eni  denote  any 

design  with  m.  observations  in  each  of  the  intervals  I,  . .  Then  E  converres 
1  !:  ,i  ni 

weakly  to  E.  Furthermore,  sup  )E  (x)  -  E(x)(  <  m./n..  So  since  lim  m./n1,/“  =  9 

X  1 1  .  11  l-Ko  j  i  ’ 

1 


12 


!/•> 

lir,  n.  (C  (x)  -  S(x))  =  0.  Dy  Corollary  4.1,  the  mean  vector  in  Theorem 
!-*•  1  ni 

4.1  is  zero  for  a  sequence  of  designs  constructed  in  this  way.  A  special 

case  of  interest  is  the  case  of  m.  =  M  for  some  fixed  N  and  n.  =  Mi. 

i  i 

Example  2 •  Cere  is  an  example  of  a  sequence  of  designs  for  which  the 
h(x)  function  is  non- zero  and,  by  Theorem  4.2,  the  mean  vector  b  in  Theorem 


4.1  is  not  zero  for  all  mean  functions  m(x).  Let  [u,  v]  =  [0,  1]  and  let 

£  be  the  uniform  distribution  on  [0,  1].  Let  =  i“  for  i  =  1,  2,...  .  Let 

■> 

L  be  the  desinn  which  tabes  one  observation  at  each  of  the  i“  -  i  points 
n . 

2 

1/i  +  k/i~,  k  =  i  -  i  -  1,  and  i  observations  at  1.  For  each  i, 


x  -  1/i  <  £  (x)  s  x  for  all  x  e  [0,  1]  so  5  (x)  converges  to  £(x).  But 

n .  n . 

1  1 


for  any  i  and  any  x  <•  (1/i,  1),  -1/i  <  5  (x)-  £(x)  <  -1/1  +  1/i''.  Thus  for  all 

x  £  (0,  1),  h(x)  =  lin.  n^-fC  (x)  -  ?(x))  =  lim  i(r  (x)  -  £(x))  =  -1.  Sunnose 
i-n»  1  n.  ^n- 

1  7  1 

f  (x)  =  1  and  f0(x)  =  x.  Sunpose  n(x)  =  x".  Then  using  the  result  of  Theorem 


4.2,  the  mean  vector  b  can  be  calculated  to  be  b  =  (-1,  2)1 
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